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Abstract 

The dyonic 1/4-BPS states in 4D string theory with Af = 4 spacetime su- 
persymmetry are counted by a Siegel modular form. The pole structure of the 
modular form leads to a contour dependence in the counting formula obscuring 
its duality invariance. We exhibit the relation between this ambiguity and the 
(dis-) appearance of bound states of 1/2-BPS configurations. Using this insight 
we propose a precise moduli-dependent contour prescription for the counting for- 
mula. We then show that the degeneracies are duality-invariant and are correctly 
adjusted at the walls of marginal stability to account for the (dis-) appearance 
of the two-centered bound states. Especially for large black holes none of these 
bound states exists at the attractor point and none of these ambiguous poles 
contributes to the counting formula. Using this fact we also propose a second, 
moduli-independent contour which counts the "immortal dyons" that are stable 
everywhere. 
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1 Introduction 

A microscopic counting formula was proposed more than a decade ago for dyonic BPS 
states in the M = 4, D = 4 theory pQ corresponding to the toroidally compactified 
heterotic string. The interest for this counting formula was revived about two years ago 
in the context of higher-order curvature corrections to the entropy [2]. Subsequently, 
the formula was given a proper derivation [3] using the 4D-5D connection [4J (See 
also [5] for a toroidal example) and the known results for the microscopic counting of 
five-dimensional black holes [6j [7]. A particular feature of the formula, namely the 
occurence of a genus two modular form, was given a novel interpretation in terms of 
string networks [8]. Recently, a similar counting formula was proposed for a class of 
more general J\f — 4, D — 4 theories [§], known as the CHL models |10j . This class of 
models has subsequently been studied from various different angles [TT] . 

In the meantime, various puzzles have been raised about these dyon counting for- 
mulas [T2l H3] . First of all, it has been observed that there is a subtlety in checking 
their S'-duality invariance. Secondly, there is an ambiguity in choosing the integration 
contour arising from the complicated pole structure of the modular forms that enter 
the formulas. Finally, it has been noted that the BPS spectrum in the macroscopic 
supergravity theory is subjected to moduli dependence due to the presence of walls of 
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marginal stability for some multi-centered bound states. See [H] and references therein 
for a discussion of this phenomenon in the Af = 2 context. Finally, either by using 
a duality argument |12j . or by studying a specific example in great details 1 1 3 J . there 
have been some hints that all the above issues might actually have something to do 
with each other. 

The goal of the present paper is to address these issues and provide a resolution to 
some of these puzzles. In particular, our aim is to present a precise contour prescription 
that will lead to a counting formula that is manifestly 5*-duality invariant. In fact, we 
will find two natural prescriptions of this kind, one moduli dependent and a second 
only depends on the charges. To arrive at these prescriptions, an important role is 
played by the one-to-one correspondence between various poles in the integrand of 
the counting formula, and the different decay channels in which a dyon can be split 
into two 1/2-BPS particles. This correspondence between poles and bound states was 
envisaged in pQ, and was recently reiterated in [TJ1 H3]. It turns out that the only 
poles that can be crossed when the choice of contour is varied are precisely the ones 
that admit such a correspondence. Moreover, we find that the contributions of the 
poles exactly match the expected number of states corresponding to the two-centered 
configurations of BPS dyon^J The key observation which allows us to identify the 
correct contour prescription is that the resulting integration contour should render the 
counting formula explicitly S'-duality invariant, and should furthermore automatically 
take the (dis-) appearance of the two-centered bound states into account when a wall of 
marginal stability is crossed. This leads to a moduli dependent degeneracy (or index-) 
formula that counts all the living dyons in every region of moduli space. In particular, 
we observe that the walls of marginal stability have the property that for large black 
hole charges (as opposed to "small black holes" with vanishing leading macroscopic 
entropy), none of the two-centered bound states of 1/2-BPS particles can exist when 
the background moduli are fixed at their attractor value. Using this fact we also 
propose a second, moduli-independent contour prescription, which has the property of 
counting only the "immortal dyons" which exist everywhere in the moduli space. 

The paper is organized as follows. In section two we start by reviewing the dyon 

3 Recently this fact was independently noted in 15], which appeared while this paper was being 
prepared. 
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counting formula and formulate the issue of its contour dependence. The contour 
dependence of the microscopic counting formula is analyzed in section three. We derive 
the condition for a given pole to contribute to the degeneracy formula and calculate its 
specific contribution. In section four we give more details of the macroscopic theory and 
derive the stability condition for the two-centered bound states of 1/2-BPS particles. 
In section five we relate the two sides and present our two contour prescriptions, one 
corresponding to the "jumping" index and one to the "eternity" index. Finally, we 
finish with some discussions and open questions in section six. 

2 The Microscopic Counting Formula and the Poles 

In this paper we are interested in dyonic BPS states in string compactifications to 
four dimensions with M = 4 space-time supersymmetry. The simplest of these string 
theories is the K3 x T 2 compactification of the type II string, or equivalently the 
toroidally compactificied heterotic string. Its [/-duality group is the product of the 
SL(2, Z) electric-magnetic duality and the 0(6, 22; Z) T-duality group. A more general 
class of four-dimensional M = 4 string theories is obtained by taking the orbifolds of 
the aforementioned theory. For definiteness, in the following we will present our results 
for the simplest case without any orbifolding, and in the end briefly comment on how, 
for a class (the CHL models [10]) of the orbifolded theories, the same steps can be 
modified and followed to arrive at a very similar result. 

The BPS states that preserve one-half of the supersymmetry are well understood. 
Using the duality symmetries, these 1/2-BPS states can be mapped to purely electrically 
charged states corresponding to heterotic strings carrying only momentum and winding 
charges. Their degeneracy follows from the level-matching condition and the state 
counting of the 24 right-moving bosonic oscillators of the heterotic string. One finds 



None of these states carry any macroscopic entropy, at least not at the leading order. 
To obtain a macroscopic entropy it is necessary to consider BPS states which preserves 




(2.1) 
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only one-quarter of the supersymmetry. Such 1/4-BPS states necessarily carry both 
electric as well as magnetic charges 



(p,Q)£r 6 , 22 ©r 6 , 22 , (2.2) 

and their leading macroscopic entropy is given b}|^] 



S = 7i\P A Q\ = n^/Q 2 P 2 - {Q ■ P) 2 , (2.3) 

where the inner product is the standard SO(6, 22)-invariant one on r 6 2 2- The counting 
formula for the 1/4-BPS states proposed in [T] takes the following form: 



p —iw ( g ) ^ ( Q ) 

D(P, Q) = j dQ (-1) (P ' Q)+1 , (2.4) 

where we have incorporated the sign factor (— l) p '^ following [19] . Here $(fi) is an au- 
tomorphic form of the genus two modular group, which means that under the Sp(2, Z) 
transformation^] 

Q -f {AQ + B){CQ + D)- 1 (2.5) 

it transforms as 

-> (det(Ctt + D)\ $(0) . (2.6) 

For the case of the toroidally compactified heterotic string the weight k is equal to 10. 
The automorphic form $(fi) is a well-defined function on the Siegel upper-half plane 
defined by 

det(Imfi) > 0, Tr(Imfi) > . (2.7) 

The precise expression for $(f2) will not be important for the purpose of our paper. The 
main property of $ that will concern us is that it has double zeroes at specific loci in 
the Siegel domain. In the entropy formula these lead to double poles in the integrand. 
When one identifies Q with the period matrix of a genus two surface, the poles in 1/$ 



4 The short hand notation \P A Q\ uses the analogy with the norm of the exterior product of two 
vectors. Note the r.h.s. indeed vanishes when P and Q are parallel. 

5 The 2x2 matrices A, B, C, D have all integer entries and satisfy the following relation 

AB T = B T A ■ CD T = DC T ; AD T - BC T = l 2x2 ■ 
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occur precisely at those values of f2 at which the genus two surface degenerates into two 
separate genus one surfaces through the pinching of a trivial homology cycle. These 
degenerations are labelled by elements of Sp(2, Z) and are characterized by the fact 
that the transformed period matrix is diagonal. We will write this condition as 

((Att + B)(Ctt + D)- 1 ^ = , (2.8) 

where the superscript fl denotes the upper right (or equivalently, lower left) element of 
the symmetric 2x2 matrix. 

The poles with C / play an important role in establishing the correspondence 
between the macroscopic and microscopic entropy. However, in this paper we will not 
be concerned with these entropy-carrying poles, instead we will restrict our attention 
mostly to the ones labelled by the elements of Sp(2, Z) with C = B = and A = 
(D T )~ 1 . These elements constitute the SX(2,Z) subgroup of Sp(2, Z) corresponding 
to the electric-magnetic or ^-duality. Under the action of the element 

7 =( a j ), ad-bc=l, (2.9) 



c d 

the charges P and Q transform into 



Pv \ f P 



7 



7 n ■ (2-10) 



Q^y J y Q 

The period matrix should transform in such a way that the exponent in the counting 
formula remains invariant, that is Q — * fi 7 with 

n 7 = ( 7 r ) -1 n 7 - 1 . (2.ii) 

Hence, the poles of the Siegel modular form $(0) corresponding to these SL(2, Z) 
elements are located at 

fi 7 = . (2.12) 

Due to the presence of these poles, one has to be careful with choosing the contour 
C: the counting formula will "jump" when the contour crosses one of these poles. 



Therefore, strictly speaking the formula (2.4) for D(P, Q) is not just a function of the 



charges P and Q but also depends on the contour. 
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One of the problems with this contour dependence is that it obscures the invariance 
of the counting formula under the S-duality. As mentioned before, the exponential 



factor in the degeneracy formula (2.4) is invariant under the simultaneaous SL(2,Z) 
transformation of the charges and Q 

e~™ (q)Hq) ^--(SMS) , (2.13) 



while we also see from the Sp(2, Z) transformation property (2.6) that the modular 
form is invariant under the transformation of its argument: Q — > f2 T Furthermore, 
using the fact P 2 ,Q 2 = mod 2 and ad+bc=l mod 2 one shows that 

(_l)( p -Q) = (_l)( p 7-Q7), (2.14) 



Therefore the integrand of the degeneracy formula (2.4) is invariant under S'-duality. 
This fact is not yet sufficient, however, to prove the invariance of the degeneracies. 
Namely, due to the presence of the poles, the expression for D(P,Q) fails to be S'- 
duality invariant, unless the contour C is also transformed to a new contour C 7 . Namely, 
the equality 



dSl -^r- ' = i ^ Km <- 1)( ™ (2 - 15) 



only holds when upon inserting (2.11 ) on the r.h.s., the new contour C 7 in the f2 7 -plane 
is the same as C in the fi-plane. 

The contours C and C 7 are really different, and in general cannot be deformed into 
one another without picking up any residue. Therefore, one concludes that under S'- 
duality the choice of contour has to change. A natural way to achieve this is to let 
the contour depend on the charges, and possibly also the moduli fields, since these 
quantities do transform under S-duality. Indeed, there is an important reason to 
suspect that the dyon counting formula is moduli-dependent, since it is known that 
certain multi-centered BPS solutions only exist in some range of background moduli 
and decay when a wall of marginal stability is crossed. The aim of this paper is to 
determine the charge and moduli dependence of the contour, so that it is consistent 
with S-duality and also takes into account the decay of dyonic bound states. 
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3 Contour Dependence and Pole Contributions 



In this section we will examine the contour dependence of the dyon counting formula 
in details. For this discussion it will be convenient to parametrize the elements of Q 
explicitly as 

V=( P . V _). (3.1) 



v a 

In terms of these quantities, the formula for the microscropic dyon degeneracies D(P, Q) 
takes the form 

r e -iTr(P 2 p+Q 2 a+(P-Q){2u+l)) 

— (b dpdodv — , (3.2) 

Jc HP,*,") 

where C is a contour in the Siegel domain defined by 

Imp > 0, Ima > 0, Implmcr > (Imz/) 2 . 

In the original proposal of p], the degeneracies were expressed in terms of the expan- 
sion coefficients of 1/$ in powers of e 2mp , e 2mcr , and e 2mu . As explained above, this 
prescription is somewhat ambiguous, since the expansion will depend on the location 
of the contour with respect to the poles. 



So let us have a closer look at the possible choice of the contour C in (3.2). Due to 
the fact that we are dealing with a modular form, the contour will have to be inside 
a fundamental domain of the Sp(2, Z) modular group. A natural choice of contour is 
to perform the integral over the real parts of p a and u, while keeping the imaginary 
parts fixed. Specifically, the integration range of the real variables is 

< Rep, Rea, Rez/ < 1 . (3.3) 

The integration contour is thus a three-torus. The location of the contour is determined 
by a choice of the imaginary parts. To make sure that 1/$ has a well-defined expansion 
up to high order, we will choose these imaginary parts so that fl lies well inside the 
Siegel upper half plane, that is 

det(Iiml) = Imp Ima - (Imu) 2 > 1 . (3.4) 
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Figure 1: (a) The Siegel upper-half plane for the modular form $ is the future light-cone in the 
Minkowski space IR 1,2 , and we consider the space of all contours to be a sheet of hyperboloid inside this 
light-cone, with all the points on the hyperboloid having a large distance from the origin, (b) A pole 
corresponding to an element 7 e S'L(2,Z) is a plane in M. ' which always intersects the hyperboloid 
along a hyperbola. 

To visualize the location of the poles relative to the contours, it is convenient to regard 
the vector 

(Imp, Imcr, Imz/) G M 1 ' 2 

as a vector in a three dimensional Minkowski space, with det(Imfi) playing the role 
of the 5*0(1,2) invariant inner product. The Siegel domain corresponds to the space 
inside the future light-cone, while the space of contours for a given large value of 
det(ImO) is identified with a sheet of a hyperboloid far out inside the future light- 
cone. This is shown in Figure (Tj. As mentioned before, all the double poles of the 
generating function 1/$ are located at divisors given by the Sp(2,Z) modular images 
of the divisor v = in the (p, a, v) space. The location of this v = divisor can 
be viewed as a hyperplane which transversely intersects the space of contours. The 
location of the other poles can be obtained by acting with the Sp(2,Z) group. Note 
that the Sp(2, Z) group action can be identified with that of the conformal group of 
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>1,2 



This fact can be used to show that these general divisors take the form 

kp + £a + mv + r (pa — u 2 ) + s = with 
k,£,m,r,s GZ , m 2 - Akt + Ars = 1. (3.5) 

The poles at divisors with r = 1 have exponentially dominant contribution to the 



degeneracy formula (3.2) compared to the rest in the case of large charges, as explained 



in the appendix of pQ. In [TJ] it was observed that the contour space (3.4) does not 
intersect any of the poles having |r| > 1. Indeed, a look at the real part of the above 
equation reveals that, since all the entries of Ref2 run between and 1, there is nothing 
to compensate the large contribution from det(ImO) 3> 1 contained in the real part 
of po — v 2 . Hence, these poles will always contribute to the degeneracy formula no 
matter which contour we choose, since they lie lower in the light-cone. Therefore, we 
never run into the danger of having a contour which crosses one of these poles. For our 
purpose of studying the contour dependence of the integral, it is therefore sufficient to 
concentrate on the poles with r = 0. 



Since we are only interested in the poles inside the real domain of integration (3.3), 
we can restrict our attention to the poles with r = s = 0. It is easily seen that these are 
the images of the pole u — under the SL{2, Z) subgroup of Sp(2, Z), and hence can 
be labelled by the group elements 7 of SL(2, Z)[j Specifically, in terms of the integral 
matrix elements a, b, c, and d of 7 we have 

k = —bd, i = —ac m = ad + be , 

where the length condition 

m 2 - Ak£ = (ad + be) 2 - Aabcd = 1 (3.6) 

follows directly from ad — be = 1. The plane inside M 1 ' 2 defined by the imaginary part 



of the equation (3.5) for this case can thus be written as 

Im(— bd p — aco + (ad + bc)u) = . (3.7) 



6 One can show this by, for example, classifying both sets of numbers (a, b, c, d) and (k, I, m) by 
their prime factorizations. 
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The length condition (3.6) implies that the normal vector to the plane is spacelike, 



and hence these planes always intersect the contour space hyperboloid (3.4) along a 
hyperbola. Therefore, each plane divides the contours into two sub-classes in a way 
analogous to the pole u—0 (see Figure [l]). Whether the corresponding poles contribute 
to the degeneracy formula for a given charge configuration will therefore depend on the 
contour we choose. 

Let us now determine the condition under which these poles contribute, and, if 
they do, calculate their contribution. We first concentrate on the double pole at v — 0. 
Near the v = divisor the generating function has the limit 

1 1111 



a, v) An 2 v 2 r\ 2 \p) r\ 2 \a) 



i + eV)). (3i 



Notice that the last two factors in the limiting expression (3.8) are exactly the generat- 



ing function for the 1/2-BPS degeneracies (2.1). By plugging the expression (3.8) into 



the degeneracy formula (3.2) and after performing the integration over the real part of 
p and a, one gets 

(_-\)PQ (■ -2-Ki{P-Q)v 

where we have made use of ( 2.1[ ). To evaluate the remaining integral over u, we first 
consider a contour with Imz/ > 0. For this case the contour is shown in the Figure [2]. 
When the charges under consideration satisfy P-Q < 0, one can deform the contour 
to the upper infinity of the cylinder (Imz/ — > oo) where the integrand is zero without 
crossing any pole. One thus concludes that the integral yields zero. On the other hand, 
in the case P-Q > 0, the contour can be moved to the lower infinity (Imz/ — > — oo) 
where the integrand is again zero, but now by doing so we pick up the contribution of 
the pole 

- 2md,(e- 2 ^ p - Q >)l =0 = -4vr 2 (P-Q), (3.9) 

where the extra minus sign comes from the fact that we are enclosing the pole in a 
clockwise direction. For the contours with Imz/ < 0, a similar argument shows that the 
pole only contributes when (P-Q) < 0, but now with the opposite sign as above due 
to the reverse orientation in which the pole is enclosed. One therefore concludes that 
the contribution of this specific pole to the degeneracy formula ( |2.4 ) is 

(-l)( p -Q)+ 1 \p.Q\d(P)d(Q) when (P-Q)lmu>0 (3.10) 
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Figure 2: In this fi gure we show how the pole located at v — contributes to the degeneracy formula 
for contours with 1m.v > 0. (a) For charges with P ■ Q < 0, one can deform the contour to the upper 
infinity of the cylinder where the integrand goes to zero without hitting the pole, (b) For charges 
with P ■ Q > 0, one can deform the contour to the lower infinity of the cylinder, and by doing so pick 
up the residue of the pole. 

and zero otherwise. The contributions of the other poles can be determined directly 
in a similar fashion. However, they are more easily obtained by making use the fact 
that they are the SL(2, Z) images of the v = pole. Together with the fact that the 



integrand is invariant under S'-duality (2.15), it follows that the double pole of 1/$ 
located at 



z/ 7 = —bdp — aca + (ad + bcjv = 



(3.11) 



gives the contribution 



vP r Q 7 +l 



Pry ■ Qy\ d{P 1 ) diQ-y) WllCU (Pry " Q IHIZ^ > 



(3.12) 



and zero otherwise. The equation (3.12) summarizes all the contour dependence in the 



degeneracy formula (3.2) 



As we will see, the jumps in the counting formula when a contour crosses one of 
the poles are related to the decay of marginally bound 1/2-BPS particles. Specifically, 



we will argue that (3.10) precisely counts the number of states associated with the 



bound state of a purely electric 1/2-BPS object and a purely magnetic 1/2-BPS object, 



while (3.12) is associated with more general dyonic bound states that are obtained 



11 



by electric-magnetic duality. This interpretation will be discussed in more details in 
section 5, after we describe the supergravity solution corresponding to these states. 



4 Dying Dyons and Walls of Marginal Stability 

The central charge in the Af = 4 supersymmetry algebra can be written as 

Z = -$=(P L ~ rQ L ) m T m ; m = 1, .., 6 , (4.1) 

V T 2 



where r = T\ + ir 2 is the usual complexified axion-dilaton field and the left-moving 
charges are given by a six-dimensional projection of the 28 dimensional charge vectors, 

P L = ^P A , Q L — /i™ Q A ; A= 1,2, ....,28. (4.2) 

Here fi™ is a 6 x 28 matrix comprised of the 6 x 22 moduli fields parametrizing the 
coset space ^)x^(22) • Here and from now on all the moduli fields are evaluated at 
spatial infinity. 

The square of the BPS mass is equal to the largest eigenvalue of Z^Z. One can 
choose the basis such that all the gamma matrices are hermitian, and one then finds 

&Z = -\P L - tQ l \ 2 1 - 2iP™Q n iXmn ■ (4.3) 

From the fact that the operator iP™Q r lY mn satisfies 

(iP?Q n L T mn ) 2 = \P L A Ql\ 2 = Ql P 2 L - (Ql ■ Pl? , (4.4) 

one concludes that Z^Z has the following two eigenvalues 

\Z P , Q \ 2 = -\Pl-tQ l \ 2 + 2\P l AQ l \ (4.5) 

and \Z' PQ \ 2 = -\P l -tQ l \ 2 -2\P L AQ L \ . (4.6) 

The complex number Zp^Q with the largest norm plays in theories with M = 4 super- 
symmetry the same role as the single central charge in M = 2 theories. In particular, 
it determines the BPS mass 

M p ,q = \Z PjQ \ 
12 



and will therefore simply be referred to as the central charge. By choosing a specific 
spinor basis one can also fix the phase of Zp } q. We denote it by ap,Q, i.e. 

Zp, q = e iap -Q \Z P , Q \ . (4.7) 

Note, however, that this phase is not unambiguously defined, since it depends on the 
choice of the spinor basis. As we will explain later, the criterion that determines 
whether two 1/2-BPS objects form a bound state can be formulated as a condition on 
the relative phase between the central charges of different objects. This relative phase 
is independent of the choice of spinor basis, even though the overall phase of the central 
charges is not. 

A related comment is the following. Consider the SX(2,Z) duality transformation 



of the charges (2.10) and the axion-dilaton moduli (4.16), this transformation has the 



effect of shifting the phase of the central charges by 

Zp, q -> e- ia ~>Z PiQ , a 7 = arg(cr + d) . (4.8) 

Again, the phase shift is independent of the charges, therefore all the relative phases 
will indeed be duality invariant. 

We are now interested in knowing when a dyonic bound state might decay. First 
we concentrate on the specific decay channel of a dyonic, 1/4-BPS state with charges 
(P, Q) splitting into two 1/2-BPS particles with charges (P, 0) and (0, Q). For this case, 
the condition for a wall of marginal stability is 

M p ,q = M Pi0 + M ,q , (4.9) 

which can be rewritten as 

\Z PtQ \ = \Z P)0 \ + \Z 0iQ \ . (4.10) 
Using the fact that the total central charge obeys 

Zp,Q = Zpfi + Z 0t q , (4-11) 

one finds that the condition of marginal stability can only be satisfied when the phases 
of the central charges are aligned. An explicit expression for this condition can be 



obtained either by determining these phases, or directly from (4.9 ) by using the explicit 
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formula for the BPS mass. Both approaches require a little bit of manipulation, and 
lead to the condition 

- + #7^ = °- ( 4 - 12 ) 

T 2 \PlAQ L \ 

The next step will be to consider the other ways in which a dyon can split into two 
1/2-BPS particles, and determine the corresponding walls of marginal stability. By 
definition a 1/2-BPS state must have degenerate eigenvalues of the operator Z^Z and 
thus have parallel electric and magnetic charges 

\Zp,q\ 2 = \Z' p , q \ 2 ^P\\Q. (4.13) 

As discussed in [T2], these 1/2-BPS decay channels can be labelled by SL(2, Z) elements 



where the two terms to the right of the equivalence sign are defined by the corresponding 
terms left of this sign. This equation shows that these bound states are basically the 
SL(2, Z) transforms of the bound state of a purely electric particle with charge P 7 and 
a magnetic particle with charge Q 7 . Using the fact that the central charge -Zp,Q is 
SL(2, Z) invariant up to a phase, one finds the condition of marginal stability to be 

\Zp,q\t — l-Zp 7 ,o + ^o,Q 7 |r 7 = |Zp 7 ,o|t 7 + |Z 0i q 7 | T7 , (4-15) 

where in the last two expressions the central charge is evaluated with the SX(2,Z) 
transformed value of the axion-dilaton fields 

CLT + b 

CT + d 

After a straightforward calculation, the position of the corresponding wall of marginal 
stability turns out to be precisely given by the SL(2, Z) image of the one for the bound 
state of the purely electric and purely magnetic 1/2-BPS states. Namely, the walls of 
marginal stability for all two-centered 1/2-BPS splits are 

— + i ( p L in L) r =°- ( 4 - ir ) 



r 7 ^ ^ . (4.16) 



7 Strictly speaking, the bound states are labeled by elements of PSX(2,Z), since states related by 
exchanging the two particle are physically equivalent. 
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As mentioned in [12], the projection of the above wall of stability from the full 134 
dimensional moduli space to the upper r 7 -plane is a straight line. But when regarded 
in the original r-plane, it is a circle for generic group elements 7. 

The meaning of the presence of a wall of marginal stability is that a BPS bound 
state of two particles exists on one side of the wall and disappears when crossing into 
the other side. After deriving the location of the walls for these bound states, we would 
like to know on which side these states are stable and on which side unstable. For this 
purpose we need more information about the corresponding supergravity solutions. 



Stability Conditions from the Supergravity Solutions 

Let us now consider the four-dimensional M = 4 supergravity theory describing the 
low energy limit of the heterotic string compactified on a six-torus. The metric part of 
a stationary solution reads 

ds 2 = -e~ 2U (dt + u ■ dx) 2 + e 2U dx 2 (4.18) 
e 2U = \V A Q\ = \/V 2 Q 2 - (V ■ Q) 2 (4.19) 



V xw = V -VQ- Q-VV , (4.20) 

where the indices are contracted using the standard SO(6, 22)-invariant 28 x 28 matrix 
tjab, for example V 2 = V a V b t]ab- 

The 56 harmonic functions appearing in the above solution are 

P. A 



i 

Q A (x)=D A + y2-^- , (4.21) 



Q 

x — X 



with the 56 constants given by the asymptotic value of 23 complex scalar fields (the 
axion-dilaton moduli r and the 22 complex moduli projected from the aforementioned 
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6 x 22 moduli) afl 

C A = -lm (e-^ 9 ^ 
dZ P . 

QpA 



D A = l^[e-^ d ^-\ , (4.22) 



where the P A 's and the Qa's denote the total charges coming from all the centers. From 
this expression one immediately sees that these coefficients satisfy Q A C A = P A D A , 
since the central charge is linear in all charges. 

For the specific two-center bound state considered earlier, the corresponding super- 
gravity solution has harmonic functions given by 

pA 



v A = c A + 



IX — Xp\ 



Qa = D a + -^— . (4.23) 

\X - Xq\ 

In this case the coordinate distance between the two centers \xp — xq\ is fixed by 
the integrability condition [16], obtained by taking the divergence of the both sides of 



(|4.20|), and reads 

F ' Q C A Q A . (4.24) 



\X P - Xq\ 

After some algebra this becomes 



P-Q \P l a Ql\ fn + P L -Q L | 9r>) 



\x P - x Q \ M PtQ \r 2 \Pl/\Ql\ 

Since the distance between the two centers is always a positive number, one finds that, 
in order for the bound state to exist, the expression on the r.h.s. must have the same 
sign as P ■ Q. Therefore the bound state only exists when 



8 By evaluating the AT = 4 central charge operator Z (4.1 1 in the eigen basis of Z*Z, one can write 
the BPS equations in a way analogous to the Af = 2 case [T5] . Only 22 complex moduli made out of 
the 6 x 22 real moduli fields play a role in the solution. It is indeed known that the Af = 4 moduli 
space locally decomposes as a product of 22 vector-, 44 hyper-, and 1 tensor- multiplet scalars in the 
Af =2 language (see, for example, [iTjl. 
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and decays when one tunes the background moduli to hit the wall where the above 
expression vanishes. More precisely, one finds that the distance between the two centers 
goes to infinity, and the bound state no longer exists as a localizable state. 

For the other bound states of 1/2-BPS particles obtained by acting with an element 
7 of the electric-magnetic duality group there exist similar solutions. But in this case 
the harmonic functions will have a seemingly more complicated form than the P\Q 
split studied above. More explicitly, now the harmonic functions and the corresponding 
integrability condition takes the form 



V A = C A + 



pA pA 

\ _ r ,,\ r \ _|_ r 2 

-£p 7 | \x-xc 

Qa,1 Qa,2 



Qa = D a + J* t i + " Z i , (4-27) 
and 

■ ^ Q \ | + dY Px - f ,^ P \ , + cY Q x = , (4.28) 
\xp y -x Q ^\ J Wxp^-xq^I J 

where P\ 7 2, Qi,2 are given in terms of the original charges and the group element 7 as 



(4.14). Plugging in the charges, and after some manipulations using (4.8), the above 
integrability can be written, as expected, in a similar form as above: 

(iVQ 7 ) \Pl A Ql\ fr 7il , (P l -Ql)j 



+ o , ;T • ( 4 - 2 9) 



Wp^-xqJ M P:Q V r 7,2 \Pl/\Ql\ 

Thus, following the same reasoning, one finds exactly the same stability condition 

T 7 ,l (P L ■ Ql)i 



but now with both the charges and the axion-dilaton transformed with 7 e SL(2, Z). 



5 The Contour Prescriptions and Their Interpretation 

Let us now return to the problem of identifying the contour that should be used in 
the counting formula, so that it counts the right number of states for a given value 
of the moduli. The key observation which will allow us to find the right prescription 
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is that the contour dependence due to the crossing of the pole labelled by 7 should 
exactly match the physical decay process of the corresponding dyonic bound state. For 
example, at the wall of marginal stability of the bound state of an electric 1/2-BPS 
particle with charge P and a magnetic 1/2-BPS particle with charge Q, one expects 
the degeneracy D(P,Q) to be adjusted by a certain amount corresponding to the 
degeneracy of this (P, 0), (0, Q) bound state. This degeneracy can be found in the 
following way [121 H31 EH] • Firstly, each of the two centers has its respective degeneracy 



d(P), d(Q), which is given by the 1/2-BPS partition function of the theory as (2.1). 
Secondly, there is an extra interaction factor due to the fact that the spacetime is no 
longer static. The conserved angular momentum, after carefully quantizing the system 
[T8] . turns out to be 

2J+1 = \P-Q\. (5.1) 

One therefore concludes that the jump in the counting formula when one crosses the 
wall of marginal stability from the stable to the unstable side is given by 

D(P, Q) -> D(P, Q) + (-lpQ) \p.Q\ d(P) d(Q) . (5.2) 



This jump in the degeneracy is precisely the contribution (3.10) that we found from the 
pole at v = 0! Similar jumps occur when one crosses the walls of marginal stability for 
the other dyonic states labelled by SL(2, Z) elements 7. These jumps are again precisely 



given by the contributions (3.12) of the poles at v 1 = 0. In terms of the contour space 
parametrized by Imp, Ima and Imu, we have shown that whether this pole contributes 
or not depends on the sign of (P 7 ■ Q 7 ) Imz/ 7 , while from the supergravity solution we 
have learned that whether this bound state exists or not depends on the sign of the 



l.h.s. of (4.30). It is therefore natural to make the identification 

7 " A W + \P L AQ L \J ' (5 ' 3) 

where A is a yet undetermined positive parameter. This equality actually constitutes 
an infinite number of equations, namely one for each element 7 G SX(2,Z). It is 
not immediately clear that all these equations can be imposed without running into 
contradictions. 
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One way to show their mutual consistency is to work out both sides of the equation, 
and observe that the left as well as the right can be written as a sum of products of 



the integers a, b, c and d as in (3.11 ). By identifying the various terms one arrives at a 



prescription that is independent of these integers. In other words, in this way one finds 



that the infinite set of equations (5.3) are equivalent to the following three conditions 

Pl-Ql 



Imz/ = —A 



T 2 



and 



Imp = A 



Imcr = A 



1 

T-2 



\PlAQl 

Ql-Ql 
PlAQl\ 

Pl-Pl 



(5.4) 
(5.5) 

r 2 + Jpf^l) • (5 - 6) 
These equations determine the location of the contour C in terms of the charges and 
moduli. To see that these equations are consistent with ^-duality invariance, we better 
use a more clever way to write them. It will turn out to be convenient to introduce 
the 2x2 matrices 

' ' 1 ' • ' (5.7) 



and 



M 



1 



1 

T"2 



-n 

J2 



Pl,Ql 



Ql-Ql 
-Pl ■ Ql 



-Pl ■ Ql 
Pl-Pl 



(5. 



\P L AQl\ 

Notice that the first matrix is given in terms of the asymptotic value of the axion-dilaton 
moduli, while the second depends on the charges and contains the asymptotic Narain 
moduli. These matrices transform in an identical fashion under the electric-magnetic 
iS-duality group, namely 



M T - (^y'Mrl- 1 



M 



Pl,Ql 



( 7 T ) X M Pl ,q l1 1 . 



(5.9) 



It is important to note that these transformation rules are the same as those of Q. We 
can now summarize the results of the previous section in terms of these matrices as 



follows. The location of the wall of marginal stability (4.17) labelled by the SL(2, Z) 
element 7 is given by the condition 



(5.10) 
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Here, the superscript fl again denotes the off-diagonal component of the 2x2 matrix 
inside the brackets. Similarly, the location of the corresponding pole z/ 7 = is given in 
terms of Q by 

(( 7 T )- 1 n 7 - 1 ) tt = o. (s.ii) 

In this way we are naturally led to the following moduli-dependent contour prescription. 
The contour is determined by specifying the value for the imaginary part of Q in terms 
of the matrices Ai T and M-p l ,q l containing the background moduli. The prescription 
reads 

lmQ = A {M T + M Pl , Ql ) (5.12) 

with A ^> 1. Here A is taken to be large to ensure that the series expansion of 1/$ 
converges rapidly. Moreover, as explained earlier, for large A the contour avoids all 



other poles except the ones given by v 1 = 0. Also note that the identification (5.12) is 
consistent with the Siegel condition, since 

det<imsi)=A2 GsyL >o < 5 - i3 > 

and the trace of Imf2 is also easily seen to have the required sign. Using the results 



(5.10) and (5.11), one easily verifies that this contour precisely crosses the right poles 
at the walls of marginal stability to account for the correct jumps in the dyon degenera- 
cies. Furthermore, note that the contour prescription leads to a manifestly ^-duality 
invariant counting formula. 



The Attractor Contour for Large Charges 

For large charges corresponding to a macroscopic black hole, it is natural to ask what 
happens to our prescription when one takes the moduli at infinity to be at the attractor 
value. Since the attractor values of the moduli are completely determined by the 
charges, this procedure leads to a degeneracy formula that is independent of the moduli. 
At the attractor point in moduli space the following equations hold for the Narain 
moduli 

PrU. = 0, QrU.=0, (5.14) 
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and the axion and dilaton are given by 

_P 1 Q _ \PAQ\ 

7"l|attr. — Q2 ' ^ I attr. — q2 ' {0.10) 

From these equations the attractor values of the matrices M. T and M.p L ,Q L are easily 
determined. One finds 

■M T |«ttr. = Mp^Q^. = M PiQ , (5.16) 

where the 2x2 matrix Mp t Q is defined by 

1 ( Q-Q -P-Q 



M, 



(5.17) 



P ' Q ~ \PAQ\ V -P-Q P'P 

Here the inner products between the charges are again defined using the moduli- 
independent SO(6, 22) invariant metric. In this way, we find that at the attractor 
point our moduli-dependent contour reduces to the following moduli-independent ex- 
pression 

lmQ = 2AM Pi Q. (5.18) 

Again the SL(2, Z) invariance is manifest, since both sides transform in the same way, 
and hence this prescription also leads to a 5-duality invariant counting formula. But 
what are the states that are being counted by this formula? In fact, we will now argue 
that these are precisely the 1/4-BPS states that are not given by the bound states of 
two 1/2-BPS particles, and therefore cannot decay. Namely, when one fixes the moduli 



to be at the attractor values, the stability condition (4.26) reduces to 



(P ■ Q)l 

—2- 7 > n 

\PAQ\ >U ' 

which can clearly never be satisfied. In other words, none of the bound states of two 
1/2-BPS particles can exist at the attractor moduli, which is a fact consistent with the 
general phenomenon that an attractor flow always flows from the stable to the unstable 
side. In this sense, our moduli-independent contour prescription leads to a counting 
formula which counts only the "immortal" dyonic states that exist everywhere in the 
moduli space. Notice further that this class of contours is not defined for charges with 
negative discriminant, since they lie outside of the Siegel domain. Furthermore, they 
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do not have an attractor point, and there is no single-centered supergravity solution 
carrying these charges. 

Finally we would like to briefly comment on the role of the number A in our pro- 



posed contours (5.12), (5.18). It can be seen as playing the role of a regulator for 



the convergence of the generating function. To see this, notice that when we take the 



contour according to our prescription (5.18), the contribution 
D(P,Q)e <ir («)"(«) 



\D(P,Q)\e 



4Att|PAQ| ^ e S e -4AS 



(5.19) 



of certain large charges to the partition function is highly suppressed when A ^> 1, and 
we are therefore left with a rapidly converging generating function. 



6 Conclusion and Discussion 

In this paper we establish the precise relation between the contour dependence of 
the microscopic formula and the presence /absence of bound states of two 1/2-BPS 
configurations in different parts of moduli space in the macroscopic supergravity theory. 
Furthermore we propose a moduli- dependent prescription for the integration contour, 
such that all these two-centered bound states are correctly counted by the counting 
formula. Therefore we arrive at the surprising and somewhat unexpected conclusion 
that the counting formula actually counts the degeneracies in all of the moduli space, 
and that there is a well-defined way to extract these degeneracies from the counting 
formula by choosing the contour appropriately. In particular, the counting formula 
has a built-in S"-duality invariance when the prescribed contour is used. Furthermore, 
for large black hole charges, we also propose a second, moduli-independent contour 
by going to the attractor value of the moduli, using which only the "immortal dyons" 
which exist everywhere in the moduli space are counted. 

But there are certainly things we do not yet understand about this counting formula. 
First of all, what are the meaning of the o ther poles which seem always to contribute? 



Poles with r > 1 in (3.5) give a correction to the index of order e s l r = e 7r ' PA ^' //r for large 
charges, which suggests that they might account for a split of charges into r pieces in 
some way. It would be nice to understand better the role of other splittings of charges 
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in the BPS spectrum. A second but not unrelated question is, what are the spacetime 
interpretation of the modular group Sp(2,Z)? Resorting to the product expression for 
the generating function [1] 

I j / j Nc(4fc£-m 2 ) 

$Tf2l ~~ e 27ri(p+o-+!/) n ( 1 — e 27ri(fcp+fcr+m^) J 
V ; {k,£,m)>0 V 7 

reveals that all the poles susceptible to contour dependence are caused by the lowest- 
lying oscillators with multiplicity c(— 1) = 2 (recall that c(n) = for n < —1), and vice 
versa. This suggests that there might be a way of rewriting the generating function 
analogous to the sum over modular images of the polar terms as in [2TH |2"T| I22J Q3j- We 
hope to put these puzzles together and return to these issues in the future. 

Another question to be asked is, how would the inclusion of higher-order corrections 
to the low energy effective action affect our macroscopic analysis? In which way does, 
if it does, the counting formula encode the information about these corrections? 

Finally, this inverse of the modular form seems to our knowledge to be the 
first example of a moduli-independent partition function, in the sense that the index 
is always summarized by the same generating function, but different expansion points 
must be used for different background moduli. This fact might shed some light on the 
enigma of the split state counting [14J, which arises due to the presence of marginal 
stability walls in a M = 2 setup. It would be very interesting to investigate whether 
some of the similar structure is also present in (some classes of) M = 2, D=4 theories. 

As for the CHL models, a dyon counting formula has been proposed for appropriate 
Zat orbifolds of the above theory for N = 2,3, 5, 7. In these theories, the rank of the 
gauge group is reduced and the S-duality group is now the following subgroup of 
SX(2,Z): 

ri(JV) = |^ ^ eSL(2,Z) ; c = mod N,a,d=l mod JvJ . (6.2) 

Moreover, the family of the contour-dependent poles of the proposed generating func- 
tion - 1 - , which is now a modular form of a subgroup of 5p(2,Z), and the ways in 
which a dyon can split into two 1/2-BPS particles, are both modified compared to the 
original theory. Nevertheless, we find that they can again both be labelled by the 
elements of the S'-duality group Ti(N), and these poles again give the same jump of 
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index as the decaying of these bound states. In particular, following the same argu- 
ments we make exactly the same proposal (5.12) for the integration contour for the 
dyon counting formula of this class of models. 
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